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The atmospheres of a large proportion of white dwarf stars are polluted by heavy elements! 
that are expected to sink out of visible layers on short timescales23, This has been inter- 
preted as a signature of ongoing accretion of debris from asteroids4, comets, and giant 
planets’, This scenario is supported by the detection of debris discs” and transits of plan- 
etary fragments®! around some white dwarfs. However, photospheric metals are only indirect 
evidence for ongoing accretion, and the inferred accretion rates and parent body composi- 
tions heavily depend on models of diffusion and mixing processes within the white dwarf 
atmosphere! Here we report a 4.40 detection of X-rays from a polluted white dwarf, 


G29-38, using a 106ks exposure with the Chandra X-ray Observatory“, demonstrating di- 


rectly that the star is currently accreting. From the measured X-ray luminosity, we find an 
instantaneous accretion rate of Mx = 1.63149 x 10° gs71. This is the first direct measure- 
ment of the accretion rate onto the white dwarf, which is independent of stellar atmosphere 
models. This rate exceeds estimates based on past studies of the photospheric abundances 
by more than a factor two, and implies that convective overshoot has to be accounted for in 
modelling the spectra of debris-accreting white dwarfs. We measure a low plasma temper- 
ature of kT = 0.5 + 0.2 keV, corroborating the predicted bombardment solution for white 
dwarfs accreting at low accretion rates3114| Offering a new method for studying evolved 
planetary systems, these observations provide the opportunity to independently measure the 


instantaneous accretion rate of planetary material, and therefore investigate the timescale of 


accretion onto white dwarfs", and the evolution and replenishment of debris disks" 


G29-38 is among the 100 closest white dwarfs!®! and has hence been subject to detailed 
studies across most wavelength ranges. The detection of an infrared excess was initially interpreted 
as arising from a brown dwarf companion, However, subsequent ultraviolet (UV) spectroscopy 
revealed trace metals in the hydrogen atmosphere of the white dwarf? which were interpreted 
as the signature of ongoing accretion from a compact dusty debris disk that formed from the tidal 
disruption of an asteroid?” Assuming that the white dwarf atmosphere is in an equilibrium between 
accretion and gravitational settling predict an accretion rate of 6.5 x 108 gs71. Some 25-50% 
of Milky Way white dwarfs are expected to be metal-polluted!! and there are now more than a 
thousand known metal-polluted stellar remnant? and yet the evidence for ongoing accretion 


remains circumstantial, based on the modelling of atmospheric abundances. 


White dwarf accretion should be accompanied with intense heating of the infalling material, 
sufficient to promote cooling via X-ray emission! This has been observed directly for white 
dwarfs accreting from stellar companions!31!4123)24) but never for a white dwarf accreting planetary 
debris. An XMM-Newton observation of G29-38 resulted in a non-detection due to the presence of 
a bright X-ray source nearby, placing an upper limi on the accretion rate of ~ 101° gs71. Similar 
upper limits from non-detections have also been derived for a handful of other metal-polluted white 
dwarfs!5) but a detection of X-rays — direct evidence of ongoing accretion — has never previously 


been achieved. 


We observed G29-38 with Chandra ACIS-S on five occasions between 22 and 27 Septem- 
ber 2020, with a total exposure time of 106.33 ks (see Extended Data Table|i] for details). The 
sky location of the X-ray photons of our observations is shown in Fig.[I} along with an image 
at optical wavelengths. The data reduction utilised three standard Chandra science bands: soft 
(0.5-1.2 keV), soft+medium (0.5-2.0 keV) and broad (0.5-7.0 keV). Within 1 arcsec of our target 
coordinates, we detected a total of five X-ray events; with four in the soft band and one in the 
medium band. Using a Bayesian approach? (see Methods), we found the 68% confidence interval 
on the source count rate in the broad band to be (2.4-6.6) x 10~° counts per second (see Extended 


Data Table] for count rates in the soft and soft+medium bands). 


Fig.[2] shows the sky location of the soft+medium band events, which are fully consistent 
with the lo uncertainty on the position from the combined Gaia EDR324 and Chandra ACIS-S 


astrometry. Considering the three science bands we computed the probability that the observed 


events could arise by chance. We used a standard aperture photometry approach (see Methods) 
with source counts of 4, 5, and 5, in the soft, soft+medium and broad bands, respectively, and 
a measured background of 0.10, 0.21 and 0.57 counts per | arcsec aperture, respectively. Using 
the Poisson distribution, we calculated the probability of receiving the observed source counts or 
higher within 1 arcsec, given the observed background count rate, and found the statistical sig- 
nificance of the source counts in the three bands to be 4.62, 4.71 and 3.640 (see Extended Data 
Table[2), indicating a high-confidence detection of X-rays from G29-38. We note that accounting 
for the astrometric uncertainty by centring the aperture on the detected source position (Fig.[2), and 
using a smaller radius of 0.5 arcsec, which is appropriate for soft sources in ACIS-S (see Methods), 


increases the detection significance to 5.65, 5.94 and 5.060 in the three science bands, respectively. 


As confirmation that the observed source counts originated from the target, and not a back- 
ground source, we utilised a source detection algorithm (see Methods) to derive a sky density of 
sources at the depth of the ACIS-S observations. In the soft, softtmedium, and broad bands, this 
analysis resulted in 7, 25, and 32 sources, respectively, setting the background source sky density 
over the whole image to psky = 3, 10 and 13 x 10~° arcsec-*. The detected sources comprise real 
astrophysical sources as well as any spurious detections arising from background counts. Interpret- 
ing this as the probability of a background source chance-aligning within 1 arcsec of our target, we 
were able to confirm the source events originated from our target with a statistical significance 
of 4.2, 3.9 and 3.80. Finally, as an independent confirmation, we also performed a Monte Carlo 
aperture photometry analysis (see Methods) within 100 arcsec of the target. In the soft band image, 


out of 100,000 test apertures, only 0.001% retrieved four events, with none retrieving more than 


four. This allows us to rule out chance alignment with a background source at 4.40. 


We conclude that the four recorded events in the standard soft science band reveal a > 40 
detection of X-rays from G29-38, with the source detected at 4.62 or 5.650, depending on the 
aperture size, and chance alignment ruled out at 4.40. The energies of the five recorded events 
within | arcsec of our target are all in the range 0.7—1.4 keV, with four below | keV. As the sen- 
sitivity of the ACIS-S detector is greatest in the range ~1-6.5 keV, and drops off steeply towards 
lower energies (see Extended Data Fig.(3), the detection of most photons below 1 keV strongly 
suggests that the X-ray emission spectrum is very soft, and thus emitted from a relatively low- 
temperature plasma. To determine the best-fitting plasma temperature, Tx, we initially adopted an 
optically-thin isothermal plasma model and tested two distinct sets of abundances: bulk Earth 
and the spectroscopically determined photospheric abundances of G29-382! The plasma temper- 
atures were consistent, with the photospheric abundances corresponding to kTx = 0.49+9:34 keV 


(see Extended Data Table[5). 


While no debris-accreting white dwarf has previously been detected at X-ray wavelengths, 
accretion onto white dwarfs at much higher rates is common in close white dwarf binaries. Our 
measured plasma temperature of 0.5 keV from a debris-accreting white dwarf is much lower than 
for white dwarfs accreting from stellar companions, which typically accrete at rates many orders 
of magnitude higher (> 10‘ gs~!) and have plasma temperatures”! in the range 5-50 keV. This is 
a robust result because the sensitivity of the ACIS-S detector increases to higher energies (peaking 


the range 1-6.5 keV) and emission with the same luminosity at higher temperatures would have 


been readily detected. Measuring a temperature an order of magnitude lower than other accreting 
white dwarfs points to a heating mechanism that is distinct from the strong stand-off shocks thought 
to heat the infalling material at higher accretion rates. A possible explanation is that the very low 
accretion rate we observe for G29-38 is not sufficient to support a stand-off shock, and the infalling 
material impacts directly onto the white dwarf surface: the previously proposed “bombardment” 
solution L4, For a white dwarf with mass Mwp = 0.6 Mo and radius Rwp = 0.0129 Ro, the 
plasma temperature arising from accretion in the bombardment scenario was predicted by ref l 
(see their equation 9) to be Tx œ~ 0.6 keV, which is consistent with the measured plasma tempera- 


ture from our observations. 


In order to determine the X-ray flux, we also adopted a more physically-motivated spectral 
model that combines emission from an optically-thin plasma at a range of temperatures (a cooling- 


flow model; see Methods for the full set of model fits). 


The best-fit model was integrated over the Chandra ACIS-S passband (0.3-7.0 keV) to de- 
rive an X-ray flux of Fx(0.3—7.0 keV) = 1.97932 x 107 ergs~!cm™?. We also perform the 
integration over a wider energy band to estimate the flux at unobserved wavelengths. From the 
Gaia EDR3 parallax, the distance to G29-38 is d = (17.53 + 0.01) pe which implies a best-fit 
X-ray luminosity of Lx (0.3—7.0keV) = 7.24778 x 10% ergs~!. This is many orders of magni- 
tude lower than observed for white dwarfs accreting from main-sequence companions (typically 


Lx = 10—103 ergs!) 


The instantaneous accretion rate can be estimated from the X-ray luminosity using 


2 
L Rwp 


ee 
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(1) 
where Rwp = 0.0129 Rọ is the white dwarf radius, Mwp = 0.6Mo is the white dwarf mass, G 
is the gravitational constant, Lx is the measured X-ray luminosity and the factor two accounts for 
50% of the emitted photons being directed back towards and absorbed by the stal The constant 
A quantifies the fraction of the total luminosity, from accretion, carried in the observed band (0.3- 
7.0 keV). To calculate the X-ray accretion rate from our observations, we take the limiting case that 
the total luminosity is equivalent to the observed X-ray luminosity (A = 1). We allow the plasma 
to be formed near the white dwarf surface, which is supported by the bombardment scenario of 
accreted material directly impacting the photosphere!4) Under these assumptions, Fig.|3|shows the 
68% confidence interval on the instantaneous accretion rate derived for the two plasma models 
(isothermal and cooling flow) and two debris abundances (bulk Earth and photospheric), across the 
range of uncertainty on the photospheric effective temperature of G29-38. Using Equation (I), the 
bulk Earth abundances and cooling flow model, we find the best-fit instantaneous accretion rate to 
be Mx(0.3—7.0 keV) = 1.631939 x 10° gs71. This is the only measurement of the instantaneous 
accretion rate of G29-38 from X-ray observations and is the first such measurement at any metal- 


polluted white dwarf. 


The accretion rate measured from the X-rays is strictly a lower limit on the true accretion rate 
owing to: (1) limited constraints on the accretion-induced flux outside the observed wavelengths, 
and (2) the possible contribution of cyclotron cooling of the post-shock plasma. These two factors 
potentially allow an accretion rate higher than the one measured from the Chandra data alone, and 
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we consider them both in turn. 


The observations provide robust constraints on hard X-ray emission (>2.0 keV) as the Chan- 
dra ACIS-S detector is most sensitive in the range 1.0-6.5 keV (see Extended Data Fig.[3), whereas 
significant emission is likely softwards of the ACIS-S bandpass (< 0.3—0.5 keV). To estimate the 
flux carried at unobserved wavelengths we integrated the best-fit cooling flow model across the 


wider energy band 0.0136—100 keV (see Methods). We found the flux in the wider integration 


band to be Fx(0.0136—100 keV) = 4.341148 x 10715 erg s7! cm~?, with errors corresponding to 
the 68% confidence interval. This represents a flux increase of a factor 2.2173 compared with 


the Chandra ACIS-S passband integration, effectively setting the parameter A = 0.57}? from 


Equation (I). 


A magnetic field on G29-38 has the potential of funneling the accretion flow towards the 
magnetic poles of the white dwarf220) and if the density in the post-shock region is sufficiently 
high, cyclotron emission at radio wavelengths could contribute to its cooling, which would imply 
that the accretion rate based on the X-ray data is underestimated. An estimate for the temperature, 
Tp, above which cyclotron emission cooling dominates was given by equation 10 from refl such 
that one can approximate the total luminosity to be Liot ~ (Tx/Tg)Lx for a thermal plasma 
with temperature Tx. The 30 upper limit on the magnetic field of G29-38 which was found to 
be 1.5kG from spectropolarimetry observations with FORS2!5) and even if the magnetic field 
of G29-38 is close to the observational detection limit, the maximum possible correction due to 


cyclotron emission cooling is predicted to be a factor 1.11? given the observed X-ray accretion 


rate (see Extended Data Fig.[7). And for magnetic fields less than ~1kG negligible cyclotron 

cooling emission is predicted. In summary, compared to the X-ray accretion rate measured in the 

ACIS-S passband, Mx = 1.631429 x 10° gs71, the true accretion rate could be higher by a factor 
+75 


1.1+9$ due to possible cyclotron radiation and a factor 2.2+7'} due to additional flux emitted at 


unobserved wavelengths (see Methods and Extended Data Table|6). 


Accretion rates at metal-polluted white dwarfs are typically inferred indirectly from spec- 
troscopic abundance measurements, coupled with white dwarf atmospheric models that quantify 
the flux of metals moving through the photospheré23!] in the accretion-diffusion scenario. Us- 
ing this approach, the time-averaged accretion rate of G29-38 has previously been inferred2!221 
to be M = (5.0 + 1.3) x 108 and (6.5 + 1.6) x 108 gs71, where the errors indicate the typi- 
cal x25% uncertainty on spectroscopic abundance measurements) Fig. [3] shows the calculated 
X-ray accretion rates, and those derived from previous spectroscopic studies. Our observations 
establish an X-ray accretion rate which agrees to within a factor 3 of those derived via the spec- 
troscopic method. The key parameters in the determination of a spectroscopic accretion rate are 
the observed spectroscopic metal abundance, the convectively-mixed mass at the surface, and dif- 
fusion timescale at the base of the mixed region (see Equation|6] in Methods). Recent results 
from 3D radiation-hydrodynamic simulations predict that the convectively-mixed mass can in- 
crease by up to 2.5 orders of magnitude when accounting for enhanced mixing due to convective 
overshootl3334| whilst the diffusion timescale can increase by up to 1.5 dex. This leads to a 
temperature-dependent increase in accretion rate inferred from spectroscopic observations due to 


convective overshoot! At the effective temperature of G29-38 (~11,500—12,000 K), convective 


overshoot is predicted to increase the accretion rate by a factor 3—4. Thus the 3D accretion rate 


is predicted to be ~ 2 x 10° gs7! 


, which is in close agreement with the derived X-ray accretion 
rate (see Fig.B). In the atmospheric layers, convective overshoot increases both the convectively- 
mixed mass and the diffusion timescale of elements settling out of the fully mixed surface layers. 
At an effective temperature of ~12,000 K the inclusion of convective overshoot increases the typi- 


cal predicted diffusion timescale from days to years, which is more consistent with the lack of time 


variability of metal lines in G29-3 %3 


Our Chandra observations of G29-38 establish white dwarfs accreting material from dis- 
integrating planetary bodies as a new class of X-ray source, demonstrating a novel observational 
technique to investigate evolved planetary systems. In particular, X-ray observations can pro- 
vide instantaneous measurements of the accretion rate onto the white dwarf, in contrast to the 
time-averaged estimates derived from model atmosphere analyses, which are, moreover, subject 
to remaining uncertainties in the assumptions about convection and diffusion. As such, repeated 
observations of the same system will enable measurements of the temporal variability of the ac- 
cretion rate from the disk onto the white dwarf, thereby providing direct observational constraints 
on the evolution of the material within the circumstellar disks24l27, and the physical mechanisms 
driving the accretion from that disk onto the white dwarf SE] The X-ray fluxes of these systems 
are low, and only a handful are within the reach of current facilities, however, large-aperture future 
X-ray missions such as the Advanced Telescope for High-ENergy Astrophysic P? (ATHENA) will 


enable systematic X-ray studies of evolved planetary systems. 
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Fig. 1| X-ray and optical imaging of G29-38. a, Image of the recorded events in the 
soft+medium (0.5—2.0 keV) band. Black, cyan, magenta and white represent 0, 1, 2, and 
3 counts, respectively. b, The PanSTARRS1 (PS1) i-band image. In both panels, the 
white square, centered on the target coordinates, has sides of length 30 arcsec. The solid 
green circle shows 1.5 arcsec radius around target coordinates, whilst the dashed green 
circle shows the target position at the PS1 epoch (J2014.89). The Chandra observations 


reveal a significant X-ray source at the expected position of G29-38. 
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Fig. 2|Sky location of X-rays from G29-38. Position of soft+medium band (0.5— 
2.0keV) X-ray events in the merged observations. The colour of each event indicates 
the measured energy of the event. Also shown are the 10 uncertainty contours on the 
target position. The innermost, red-dashed ellipse gives the target positional uncertainty 
from Gaia EDR3. The wider, grey-dashed circle indicates the 68% positional uncertainty 
of 0.5 arcsec, arising from the astrometric accuracy of Chandra ACIS-S (see figure 5.4 of 
the Chandra Proposers’ Observatory Guide Version 23.0"). The outer, red-solid ellipse 
shows the combined 68% uncertainty interval on the position. The five source photons 


used in our statistical analyses are all consistent with the expected sky location of G29- 


38. 
“https://cxc.harvard.edu/proposer/POG/pdf/MPOG. pdf 
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Fig. 3| Accretion rates inferred from the measured X-ray flux. Accretion rates from 
the isothermal and cooling flow models are shown as open circles and diamonds, re- 
spectively. We show results for two plasma abundances; bulk Earth?” and photospherid2" 
shown with blue and red, respectively. The filled horizontal bands show the 68% confi- 
dence interval on the four accretion rates. The X-ray accretion rates are computed using 
Equation (1). Using Equation (6), the accretion rates inferred from spectroscopic observa- 
tions are shown in solid lines for a 0.6. white dwarf, assuming bulk Earth composition; 
such that the observed calcium abundance? of log (Ca/H] = —6.58 + 0.12 makes up 1.6% 
of the accreted material. Accretion rates from models which include or omit convective 
overshoot mixing” are shown in green or magenta, respectively, with 10 uncertainties 
(dashed lines) propagated from the spectroscopic abundance measurement. Also shown 
in solid circles are the previously published inferred accretion rates for G29-38 in blue=2! 


and orange?” 
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Methods 


Observations and data reduction. Our observations consisted of five exposures using Chandra 
ACIS-S carried out between 22 and 27 September 2020. Each observation had an exposure time 
between ~15-25 ks (see Extended Data Table|I), with a total, merged exposure time of 106.33 ks, 
or 29.5hr. All observations were carried out with the ACIS-S instrument using the Gaia DR222! 
proper motion-corrected coordinates of our target positioned on the S3 chip. The predicted X-ray 
luminosity was sufficiently low that our setup included the very faint mode which stores event 
grades in regions of 5 x 5 pixels, rather than the standard 3 x 3. This allows an improved identifi- 


cation and rejection of background events, particularly at hard and very soft energies. 


The data reduction was performed using the software package Chandra Interactive Analysis 


of Observations (CIAC%, The reduction began by reprocessing the five observations using the 
chandra_repro package with standard grade, status and good time filters and vfaint back- 
ground cleaning applied. We employed merge_obs to merge the five observations into a single 
set of event files along with images, point spread function (PSF) maps and exposure maps. We 
performed this routine using combinations of the standard Chandra science energy bands: soft 
(0.5-1.2 keV), soft+medium (0.5-2.0 keV) and broad (0.5-7.0 keV). During each merging rou- 
tine, the PSF maps and exposure maps were generated for each band with exposure evaluated on 


0.9 keV, which represents an approximate mean of the suspected source photon energies. 


Statistical significance of detection To determine the statistical significance of the detection, 


we utilised the Poisson distribution, as appropriate for counting statistics with small numbers of 
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events“! This allows us to quantify the confidence with which we can rule out the null hypothesis 
that we detected no source photons. We defined a source region of | arcsec around the coordinates 
of our target — Right Ascension and Declination of (352.196185 + 0.00004)° and (+05.24686 + 
0.00003)°, respectively — which have been corrected to the J2020.73 epoch using the Gaia EDR324 
astrometry. The nearby background region was defined as a circle with a 52 arcsec radius, avoiding 
obvious point sources — such as the nearby bright source which contaminated previous XMM- 
Newton observations“! (e.g., see Extended Data Fig.|I| panels a—c) — as well as drops in sensitivity 


in the exposure map. 


Extended Data Table/2| shows the source and background counts from the full, merged obser- 
vation of 106.33 ks. The Poisson distribution gives the probability of receiving N or more counts, 
given an expected mean background, b, as 


bN e~’ 


P(N|b) = = 


(2) 


We can therefore reject the null hypothesis that there are no source photons with a confidence 
level (CL) of CL = 100 x (1 — P(N|b)). In the soft, soft+medium and broad bands, we find 
CL =99.99962%, 99.99975%, and 99.97255%, respectively. The statistical significance, s, can 
be computed as 


s =erf-! (1 — P(N|b)) v2, (3) 


which for each of the three bands is s = 4.62, 4.71, and 3.640, respectively. We note that adopting 
a smaller aperture with a radius of 0.5 arcsec, which is appropriate for soft sources with ACIS-S 


(see figure 6.10 of the Chandra Proposers’ Observatory Guide Version 23.7] where the enclosed 


Thttps://cxc.harvard.edu/proposer/POG/pdf/MPOG. pdf 


20 


fractional power is >50% for 0.5 arcsec apertures), and accounting for the astrometric uncertainty 
by centring the aperture on the locus of photons at the source position (see Fig.{2), increases the 
detection significance in the three science bands to 5.65, 5.94 and 5.060, respectively. We conclude 
that in all three science bands considered, an X-ray source is detected at the expected position of 


our target, G29-38. 


Count rates and confidence intervals For each science band considered, we converted the num- 
ber of counts to a count rate using the total exposure time of 106.33 ks. In the soft, softtmedium 
and broad bands we find count rates of 0.038+9:917, 0.04710 053 and 0.04719:92f ks! per 1 arcsec 
aperture, respectively. The 68% confidence interval on count rate across all three bands is (2.1— 
7.0) x 10~° counts per second, while the 90% confidence interval spans (1.3-9.0) x 10~° counts per 
second. The confidence intervals were calculated using a Bayesian approach to Poisson statistics, 


following methodology in refl This uses a simple prior which requires a non-negative number 


of source counts. 


Investigating positional uncertainty The five recorded soft+medium band events within 1 arcsec 
of our target coordinates have a sky position which is shifted by approximately —0.5 arcsec in right 
ascension. In Fig.[2|we plot the sky position of the recorded events in the soft+medium band within 
1.1 arcsec of our target. We find that all five of the photons in the soft+medium bands lie within lo 
of the target’s position on the CCD. This confidence interval comprises two sources of uncertainty. 
The first is from the four Gaia astrometric parameters, i.e., coordinates (R.A. and DEC) and proper 


motion (Hra and Haec). The positional error at the epoch of our Chandra observations, x, is thus 
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propagated as 


OX; = (tõu) + (529 4)"]/? 3 (4) 


where t is the elapsed time between the Gaia epoch and the Chandra observations, u is the proper 


motion, Xp is the Gaia-measured coordinate, and the subscript 7 indicates either R.A. or DEC. 


The second positional uncertainty arises from the astrometric uncertainty of the telescope 
itself. From figure 5.4 of the Chandra Proposers’ Observatory Guide (Version 23.0}'| the 68% 
confidence limit on the radial offset of a given target is 0.5 arcsec. This value was computed for 
observations between 2015-2020, by comparing the radial offset of Chandra sources with optical 
sources predominantly from the Tycho-2 catalog. Combining these two sources of positional un- 
certainty, the thick ellipse in Fig.[2| shows the 68% confidence interval on the position of our target 
and illustrates the expected precision of the target position in our data. We find that all five events 


considered in our earlier statistics fall within this area. 


Ruling out background contamination To determine the probability of chance alignment with a 
background source we follow two methodologies. In the first, we run a source-detection algorithm 


to provide an independent constraint on the number of sources at this depth of pointing. We rely 


upon wavdetect#) a wavelet-based algorithm included in the CIAO package, designed for the 
spatial analysis of Poisson count data. As recommended in the documentation, we set the signifi- 
cance threshold parameter, sigthresh, to 1/npx ~ 1076, which implies that approximately one 
identified source will be a false detection. We perform this analysis in all three of the standard 


bands considered in this study, as well was for two values of sigthresh, sı and sg. In the first, 


*https://cxc.harvard.edu/proposer/POG/pdf/MPOG. pdf 
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we use sı = 10~°, which is the recommended value for the number of pixels (1024x 1024) of 
the S3 CCD. The second, more stringent, value of sigthresh, s2 = 5 x 1077, accounts for 
the ‘dead’ corners of the image which arise from aligning the CCD image with the world coor- 
dinate system (WCS). This results in an image with 1414x1401 pixels. The results are given in 
Extended Data Table|3| and Extended Data Fig.[I|shows the sky location of the detected sources in 
all three science bands. Across the entire CCD image we detect 7, 25 and 32 sources, in the soft, 
soft+medium and broad bands, respectively. We find our target is not detected in the broad band 
image, regardless of the sigthresh value. In the soft band, the target is detected at the higher 
sigthresh value, but not with the more stringent value. And in the soft+medium band we find 
the source is detected in both cases. Dividing the total number of sources in each band by the field 
of view (1024? x 0.4920? arcsec”) reveals sky densities of 2.76, 9.85 and 12.6 x 10~° arcsec” ?, 
in the soft, softtmedium and broad bands, respectively. Interpreting this as the probability of 
chance alignment with a background source, we can rule out chance alignment at a significance of 


erft (1 — Pesky) V2 = 4.2, 3.9 and 3.80, or confidence of 99.997, 99.990 and 99.987%, respectively. 


To further test this result, we performed a Monte Carlo aperture photometry experiment. 
We placed 100,000 apertures with | arcsec radii within 100 arcsec of the target. This area was 
chosen to confine the study to the region of the ACIS-S3 CCD where the PSF does not exhibit a 
significant gradient, as beyond ~120 arcsec the PSF increases significantly. This effect can be seen 
in the output of the wavdetect algorithm in the top panels of Extended Data Fig.{I] where the 


sources detected near the edge of the CCD appear to be much larger. For the Monte Carlo aperture 


photometry test we make use of the Pyt hon package sc ipy“ in particular the KDTree module. 
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This facilitates the rapid comparison of a large number of test aperture positions with those of 
recorded events from our observation. To estimate the true sky density of background sources 
we exclude the target and known nearby bright source using circular masks of radius | arcsec and 
2.5 arcsec, respectively. The final number of test apertures was slightly reduced after the removal 


of the 150-200 which fell within masked regions. 


The results of this test on the soft and softtmedium band events are shown in Extended 
Data Fig.[2| In the soft band analysis, we find 1 in 99,828 test apertures (0.001%) retrieved four 
counts, while no test apertures found more than four counts. This allows us to rule out chance 
alignment with a background source at a significance of 4.40. In the soft+medium band, we find 
56 in 99,852 test apertures (0.056%) have five or more counts, while 48 (0.048%) returned between 
6-27 counts, allowing us to rule out chance alignment at the 3.40 level. The reason the significance 
is lower in the soft+medium band is not because there is any doubt about the source detection, but 
because (1) our source has a soft spectrum and so contributes little flux in the medium band, and 
(2) there are a much larger number of background sources in the medium band, since most X-ray 
sources are harder than our target, and because the effective area of the telescope is higher. This 
test also provides an empirical confirmation of the quoted expected backgrounds in Extended Data 
Table[2} with fewer than 10% and 20% of 1 arcsec test apertures returning any counts in the soft 


and soft+medium bands, respectively. 


We have shown that the counts detected within | arcsec of our target coordinates allow the 
rejection of the null hypothesis — that no source counts were measured — with a conservative confi- 


dence of 99.97255—99.99962% (3.64-4.71c). Using the recommended aperture size of 0.5 arcsec, 
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and centring the aperture on the observed source position, we find the significance increases to 
5.06-5.910. We have also shown that the observed position of the five recorded events is consis- 
tent with the expected position of our target, to within the lo astrometric uncertainty. We have 
tested the hypothesis that recorded source counts could have originated instead from a background 
source, and by constraining the sky density of sources at this depth of pointing using wavdetect, 
we rule out a chance alignment with a background source at a confidence of 99.997 and 99.990% 
(4.2 and 3.90) in the soft and soft+tmedium bands, respectively. Confirming these results with 
Monte Carlo aperture photometry, in the soft band, with 0.001% of test apertures returning four 
counts, chance alignment can be ruled out with a confidence of 99.999% (4.40). Given the high 
degree of confidence to which the observed events can be attributed to our target, G29—38, in the 


following we perform spectral modelling to derive an X-ray flux, luminosity and accretion rate. 


Spectral model and debris composition We perform a spectral analysis assuming two differ- 
ent optically-thin plasma models, and three distinct abundance profiles. The first model is a 
one-component, isothermal plasma implemented in the vvapec model within the XSPEC soft- 
ware package (version 12.11.1)#! which uses the AtomDB atomic database“. The second model, 
mkcflow, is a cooling flow model that allows for a range of temperatures, with the relative emis- 
sion measure for each temperature weighted by the inverse of its emissivity. We selected the op- 
tion to use the AtomDB database. This cooling flow model is more physically-motivated than the 
isothermal model and is likely to provide a better estimate of the X-ray flux beyond the observed 
bandpass. We fitted the mkc £ Low model by fixing the lowest temperature plasma at the lower limit 


of the model of kT = 0.08 keV, only allowing the upper temperature of the range to vary, enabling 
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a one-parameter fit. This models the emission from material heated to the upper temperature and 
then cooling to temperatures lower than can be detected with Chandra ACIS-S. The three abun- 
dance profiles we use are Sola”, bulk Earth?” and the observed photospheric metal abundances2!! 
of G29-38 with an equal number abundance of hydrogen. The composition of the infalling mate- 
rial is best described as a rocky, water-depleted, chondritic object24) with detected lithophile (O, Si, 
Mg, Ca, Ti, and Cr), siderophile (Fe), and atmophile (C) elements. If the heated plasma is formed 
sufficiently close to the stellar surface, such as in the bombardment scenario, the rocky accreted 
material may be mixed with photospheric hydrogen. However, we found that our results were not 
sensitive to the hydrogen abundance, as at the best-fit plasma temperatures (0.5 keV) the cooling 
is dominated by metal line emission. At the distance to G29-38, the interstellar column density is 
expected“*l to be only around Ny = 5.4 x 10'8 cm~2, which has a negligible effect in the Chandra 
ACIS-S passband (<0.5% absorption at 0.5 keV). The five events in the ACIS-S spectrum were 


fitted unbinned, using the C-statistic “2 and with no background subtraction. 


Extended Data Figure |4] shows the spectral fits and Extended Data Table|5] shows the best-fit 
plasma temperatures for all six models, as well as the 68% and 90% confidence intervals. The 
bulk Earth and photospheric models agree to within the 68% confidence interval. The observed 
photospheric abundances of G29-38 could be scaled by the microscopic diffusion timescale in the 
atmosphere to infer a more accurate accreted debris composition. The difference, however, would 
be small, and the agreement between the bulk Earth and photospheric models strongly suggests 


that the observations would not be sensitive to such a correction. 
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Deriving X-ray Flux We derive a total flux due to accretion by integrating the best-fit spectral 
models over a finite frequency (or energy) range. We find a robust lower limit on the X-ray flux 
of Fx(0.3—7.0 keV) = 1.971155 x 10715 ergs~! cm~? by performing the integration only over 
the energies within the Chandra ACIS-S passband (0.3-7.0 keV). We also perform the integra- 
tion across a slightly narrower band (0.5-7.0 keV) as the instrument sensitivity below 0.5 keV has 
degraded since launch and is now relatively low. These results, with their associated confidence 


intervals, are shown in Extended Data Table|6] 


With the lower limit on X-ray flux tightly constrained, providing an upper limit on the X-ray 
flux is more challenging. This is primarily due to the lack of observations between the very soft X- 
rays (~0.1 keV) and the UV. There are no instruments currently equipped to perform observations 
at these Extreme Ultraviolet (EUV) wavelengths, so directly measuring the flux in this regime is, 
for now, impossible. The constant A included in Equation (1) describes the predicted fraction of 
the total luminosity carried in the Chandra ACIS-S passband. Previous studied41| have made the 
approximation that A = 0.5, or 0.25. We provide an upper estimate of the X-ray and EUV flux by 
integrating the best-fit spectral models over a much wider energy range (0.0136—100.0 keV). Model 
spectra over this broadband are plotted in Extended Data Fig.[5| and we include fluxes calculated 


for this broad band in Extended Data Table|6} 


Extended Data Fig.|6| shows the X-ray flux computed in the three integration bands used in 
this analysis. For clarity, from here onwards we show results only from the cooling flow model 
(mkc flow) using the photospheric abundances. The use of the cooling flow model is designed to 


provide a realistic temperature distribution and hence realistic fluxes in the EUV band (it can be 
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compared with the isothermal model in Extended Data Fig.[5p. From left to right in Extended Data 
Fig.[6]the size of the integration domain increases, and as expected the best-fit X-ray flux (and thus 
accretion rate) also increases, with the widest band (0.0136—100 keV) providing an upper estimate 
of the X-ray and EUV flux. The wider integration band results in a best-fit accretion rate which 
is 2.2*('3 times higher than that derived from the ACIS-S passband, where the errors represent the 
68% confidence interval. We note that the 68% and 90% upper bounds on the fluxes increase by a 


factor ~10 and ~100, respectively, due to the uncertain importance of the EUV band. 


X-ray luminosity We derive an X-ray luminosity from the X-ray flux measured in the previous 
section using Ly = 47d? Fx with the distance to G29-38 of d = (17.53 + 0.01) pe calculated from 
the Gaia EDR3 parallax. Extended Data Table[7] shows the best-fit X-ray luminosity for each of the 
energy ranges. From the robust lower limit on the X-ray flux, the best-fit X-ray luminosity from 
our observations is found to be Ly(0.3 — 7.0 keV) = 7.24*7-8 x 10? erg s71, where the errors are 


indicative of the 68% confidence interval. 


X-ray accretion rates Here we estimate the instantaneous accretion rate using a simple model 
which has been employed in studies of accreting white dwarfs!524#| In this model, infalling 
material reaches near-free fall velocity, forming a heated plasma as the material approaches the 
white dwarf surfac?! The plasma must radiate outwards a total energy equivalent to the initial 


gravitational potential energy of the accreted material such that24l 


Liot = 5—p > (5) 
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where Rwp is the white dwarf radius, at which infalling material impacts the atmosphere, and 
the factor 1/2 accounts for 50% of emitted photons being directed back towards and absorbed 
by the stak] We allow the emitting plasma to be formed by individual atoms reaching the 
white dwarf atmosphere in a scenario termed the “bombardment” solution, relevant for low ac- 
cretion rate systems EEA such that the X-ray emitting plasma is formed at the white dwarf radius 
(Rwp = 0.0129R.). This scenario was hypothesised for white dwarfs accreting from main- 
sequence companions at low accretion rates, and predicts a plasma temperature (~ 0.6 keV) com- 
parable to that which we have measured from our observations (0.5 keV). From the total luminosity, 
the X-ray luminosity can be written Ly = AL;,,, where the constant A accounts for the fraction 
of the flux emitted outside the observed passband. If the plasma cooling is entirely mediated by 
line cooling in the observed X-ray passband (0.3-7.0 keV), then A = 1. If the plasma experiences 
additional line cooling at harder, or softer energies than those observed, then A < 1. If the plasma 
cools via other physical mechanisms, such as cyclotron emission cooling, then some of the total 
luminosity will be radiated at radio wavelengths, in which case A < 1. This has typically been 
set to A = 0.25—0.5 in previous X-ray studies. The flux integration over the wider energy 
band (0.0136—-100 keV) allowed us to constrain the increase due to unobserved flux to a factor 
2.273, which corresponds to A = 0.51??, consistent with estimates from previous studies for 
white dwarfs accreting from main-sequence companions, with errors corresponding to the 68% 


confidence interval. 


The accretion rate can be inferred from the Chandra ACIS-S observations using Equation (I), 


which describes an accretion flow converting its gravitational potential energy into an X-ray lu- 
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minosity. In Extended Data Table|8] we use this equation in the limiting case, with all X-ray 
flux emitted within the Chandra ACIS-S passband (A = 1) to transform the X-ray luminosi- 
ties into accretion rates. We find the X-ray accretion rate, measured from our observations, to be 
Mx(0.3-7.0 keV) = 1.63*3:38 x 10° gs71. The true accretion rate could be inferred to be higher 
for two reasons: firstly, there may be additional flux carried at lower, unobserved energies, and 
secondly, there could be some contribution if the magnetic field is sufficiently high to promote 
some cooling via cyclotron emission cooling!5) We have constrained the first consideration by in- 
tegrating the best-fit model spectra over the full EUV and X-ray regime, finding the flux to increase 
by up to a factor 2.2175, or A = 0.51). We have robust constraints on the second consideration, 
finding a maximum increase due to cyclotron emission cooling of a factor 1.1196 (see Extended 
Data Fig.{7). Combining both possibilities for an increase in accretion rate, the upper estimate on 


the true accretion rate could be a factor 2r higher than the observed X-ray accretion rate, where 


the errors represent the 68% confidence interval. 


Atmospheric parameters G29-38 was identified as a metal-polluted, hydrogen atmosphere (DAZ) 
white dwarf! It is also a pulsating, or ZZ Ceti, star with large-amplitude non-radial pulsa- 
tions with periods on the time scale t ~ 110 — 1250 SOS) From fitting the broadband energy 
distribution®=23) the effective temperature of G29-38 has been estimated to be Tag = (11529 + 
206) K and (12090 + 229) K, when using Gaia EDR3 and PanSTARRS1 (PS1) photometry, re- 
spectively. From the Gaia EDR3 parallax, the mass was found to be M = (0.629+0.016) Mo and 
(0.670 0.017) Mo, when using Gaia EDR3 and PS1 photometry, respectively. This spans the full 


range of effective temperatures previously derived for this star using the spectroscopic method to 
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fit the Balmer absorption lined PEEB Using the observed Balmer line spectrum time-averaged 
over multiple pulsation cycles and the latest model spectra)! with 3D corrections= we obtain 
Teg = (11906 + 190) K and M = (0.689 + 0.031) Mo. The model accretion rates (solid lines) de- 
rived in Fig.B]are plotted over an x-axis range to include both the EDR3 and PS1 effective temper- 
atures, for a H-atmosphere white dwarf with a mass of M = 0.60 Mo. Radiation-hydrodynamic 
simulations of convective overshoot are currently only available for H-atmosphere white dwarfs 
with surface gravities of logg = 8.0 (Mwp = 0.6Mọ). For self-consistency, the calculation of 
X-ray accretion rates (Equation{I) and spectroscopic accretion rates (Equation|6) use the canonical 
white dwarf mass of Mwp = 0.6Mọ. This is 3—15% smaller than the mass of G29-38 derived 
from photometric and spectroscopic fits (Wwp = 0.629—0.670M.). Given that the X-ray accre- 
tion rate scales proportionally with white dwarf mass, this represents only a small uncertainty in 


the measured accretion rate. 


Spectroscopic accretion rates In the accretion-diffusion scenario, under the steady-state assump- 
tion, an accretion rate can be inferred from spectroscopic observations abi 


p Mevz 
Mi = Xi , (6) 


taitt,i 


where X; is the photospheric abundance of element 7, Mevz is the convectively-mixed mass and 
tain, is the diffusion timescale set by the microscopic physics at the base of the mixed surface 
layers. From spectroscopic observations and model atmosphere analyses, the time-averaged ac- 
cretion rate of G29-38 has previously been inferred122) to be M = (5.0 + 1.3) x 108 and 
(6.5 + 1.6) x 10°gs~'. The uncertainty on the first valu? was given as a typical 25% uncer- 
tainty (0.1 dex) on the spectroscopically measured calcium abundance. The error on the second 
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valu" was also estimated in the same way. Our observations establish an X-ray accretion rate 
which agrees to within a factor 3 of previously derived spectroscopic accretion rates. Recent re- 
sults from 3D radiation-hydrodynamic simulations predict a temperature-dependent increase in 
accretion rate inferred from spectroscopic observations due to convective overshoot At the ef- 
fective temperature of G29-38 (~11,500-12,000 K), the accretion rate correction is predicted to 
increase by a factor 3-4. Thus the 3D accretion rate is be predicted to be ~ 2 x 10° gs~!, which 


is in agreement with the derived X-ray accretion rate (see Fig.3). 
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Data Availability The data that support the plots within this paper and other findings of this study are 


available from the Chandra Data Archive. The observation ID numbers are given in Extended Data Tablefi] 


Code Availability The official Chandra reduction software package CIAO — which includes merge_obs, 


wavdetect and XPSEC — is freely and publicly available (cxc.cfa.harvard.edu/ciao/). So too 


is the Python package scipy. 
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Extended Data Fig. 1|Estimate of sky density from a source detection algorithm. 


Output of the wavdetect source detection algorithm using the recommended signifi- 
cance threshold, sigthresh, of sı = 1076. From left-to-right are the results for the 
science bands used in this study; soft, softt+medium and broad, with 1, 2, and 3 counts 
shown in cyan, magenta, and white, respectively. The top panels show the full field of view 
of the S3 chip on ACIS-S with sources identified by wavdetect shown in green. The bot- 
tom panels give a magnified view of the vicinity near the target, where the white square 


has sides of length 30 arcsec and is centered on the target coordinates. The source at the 
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sky location of G29-38 is detected in the soft and soft+medium band images. The number 


of sources and corresponding sky density for each band can be found in Extended Data 


Table[3] 
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Extended Data Fig. 2|Monte Carlo aperture photometry. a & c, The blue points 
show the positions of the ~100,000 test apertures, each 1 arcsec in radius, used to sam- 


ple 52 arcsec around the target. The absolute number of test apertures, after removing 
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those that fell within a masked region, is shown in the panels. The sky coordinates of 
all recorded events that fall within a masked region are shown in orange. b & d, The 
normalised histogram shows the fraction of test apertures with event counts equal to or 
greater than that of a given bin. The Monte Carlo was performed on the soft (a-b) and 
soft+medium (c—d) bands. The soft band analysis has 0.001% of test apertures returning 


four counts, allowing us to rule out chance alignment at 4.40. 
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Extended Data Fig. 3 | Spectral modelling of observed X-ray events. The observed 
counts are shown along with the best-fit spectral model for the isothermal plasma with 
photospheric abundances"! (dark blue). The effective area of the ACIS-S detector is 
shown in green, hatch. The absence of harder X-ray events (>2.0 keV) in the Chandra 


observations demonstrates that the plasma emission spectrum is very soft. 
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Extended Data Fig. 4| Best-fit spectra for the observations using XSPEC. Bottom, 
in units of instrumental counts we show the five recorded events (black) and six best-fit 
spectral models assuming Solar (magenta), bulk Earth?” (blue) and photospheri¢2#! (red) 
abundances, with either the vvapec isothermal (solid) and mkc flow cooling flow (dotted) 
plasma models. We also indicate the dominant metal emission lines (O, Mg, Si and Fe) 
from the isothermal, photospheric abundance model. Top, in real flux units, we show the 
synthetic spectra for the photospheric abundances with the isothermal (solid) and cooling 
flow (dotted) plasma models. The modelling suggests the most likely origin of the source 


photon at 1.3 keV was a Mg transition. 
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Extended Data Fig. 5| Estimate for flux carried at XUV wavelengths. Spectral en- 
ergy distribution of the best-fit isothermal (blue) and cooling flow (orange) plasma models 
with bulk Earth abundances?” down to the extreme ultraviolet (EUV) energy regime. Also 
shown are the standard Chandra science bands; soft, medium and hard. Both models 
provide a convergent fit within the Chandra ACIS-S passband, but the cooling flow pro- 


vides a more physical and larger estimate of the lower-energy flux. 
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Extended Data Fig. 6 | X-ray flux and accretion rate. X-ray flux measured in 3 bands: 
(a) 0.5-7.0 keV, (b) 0.3-7.0 keV and (c) 0.0136—100 keV, using the cooling flow model for 
the photospheric abundances?" are shown in open diamonds. The filled horizontal bands 
show the 68% and 90% confidence intervals on the X-ray accretion rate, with bounds 
shown in dashed and dotted lines, respectively. The X-ray accretion rates are computed 
using Equation (i), with A = 1, Rwp = 0.0129Ro, and Mwp = 0.6Mọ. The accretion 
rates inferred from a calcium abundance of log [Ca/H] = —6.58 + 0.12, and assuming 
bulk Earth composition (such that calcium accounts of 1.6% of the accreted material), 
are shown in solid lines. The models including and omitting convective overshoot mixing 
are shown in green and pink, respectively. Also shown in solid circles (blue and orange) 
are the previously published inferred accretion rates for G29-38, based on photospheric 


abundances from spectroscopic observations2#l221 
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Extended Data Fig. 7 | Limit on cyclotron emission cooling as source of additional 
luminosity. a, An estimate for the total luminosity from the measured X-ray luminosity, ac- 
counting for cyclotron emission cooling. We compare the measured plasma temperature, 
kTx = (0.61 + 0.28) keV, given by the cooling flow model and photospheric abundances, 
with the critical plasma temperature, Tg, above which cyclotron emission cooling domi- 
nates, defined by equation 10 from ref, The authors provide the ratio Tx/Tg ~ Liot/ Lx 
as an estimate of, for a range of accretion rates and global magnetic field strengths, the 
predicted increase in total luminosity compared to X-ray luminosity if the plasma tempera- 
ture is sufficient to be dominated by cyclotron emission cooling. The horizontal dotted line 
indicates Tx/Tg = 1, where no correction is expected below this. The vertical dotted line 
indicates the 3c upper limit on the magnetic field strength from FORS2 spectropolarimet- 
ric observations! The solid lines indicate the increase in total luminosity when compared 
to the observed X-ray luminosity. b, Predicted additional luminosity for an assumed global 
magnetic field at the 3c limit (1.5kG) across the full range of plasma temperatures and 


accretion rates calculated in this work (see Extended Data Tables|5|&(8). White space 
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indicates no additional luminosity. The upper plasma temperature from the cooling flow 
model and accretion rate derived from the isothermal plasma model is shown (solid) along 
with the 68% uncertainty (dashed). Even at the observational upper limit, the predicted 


increase due to cyclotron emission cooling is a factor of 1.1709. 
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Extended Data Table 1 | Chandra observations of G29-38. 


Obs ID Instrument Exposure Time Start Date Epoch 
(ks) (yyyy-mm-dd) (yr) 
24257 = ACIS-S 24.58 2020-09-22 J2020.727 
24256 ACIS-S 14.89 2020-09-24 J2020.732 
24658 = ACIS-S 14.89 2020-09-24 J2020.732 
23379  ACIS-S 26.23 2020-09-26 J2020.738 
24657 = ACIS-S 25.74 2020-09-27 J2020.740 
Total: = 106.33 = = 


Details of Chandra observations carried out between 22 and 27 September 2020 with 


G29-28 as the target (Pl: Cunningham). 
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Extended Data Table 2| Statistical significance of source detection. 


Band Energy Source Background CL Significance 
(keV) (per ap.) (per ap.) (%) (c) 
T1.0 To.5 T1.0 To.5 T1.0 T0.5 
Soft 0.5-1.2 4 0.100 0.025 99.99962 99.9999984 4.62 5.65 
Soft+Medium 0.5-2.0 5 0.205 0.051 99.99975 99.9999997 4.71 5.94 
Broad 0.5-7.0 5 0.566 0.142 99.97255 99.9999590 3.64 5.06 


Source counts (N) and expected background (b) in the three standard ACIS energy bands 
used in this study, given in units of counts per 1 arcsec aperture (r1.o). From top to bottom, 
the total background counts in each band was N,=270, 555 and 1530 (in the 52 arcsec 


radius background region). The expected background in the 1 arcsec radius source region 


was thus computed as b = N;,/52?. This includes reprocessing using the VFAINT mode 
for background cleaning. For an expected background, b, the Poisson distribution gives 
the probability of receiving N or more counts as P(N |b) = bNe~®/(.N!). We therefore reject 
the null hypothesis of detecting no source photons with a confidence of CL = 100 x (1 — 
P(N|b)). We also provide results for the relocalised 0.5 arcsec aperture (70.5), and note 


that the source counts are the same for both aperture sizes. 
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Extended Data Table 3| Sky density of sources from the wavdetect source detec- 


tion algorithm. 


Band No. Sources Target? Sky Density CLaign Significance 
(10-° arcsec?) (%) (o) 
Sı S2 S1 S2 Sı S2 Sı Sı 
Soft 7 6 vx 2./6 2.36 99.99724 4.2 
Soft+Medium 25 23 V xy 9.85 9.06 99.99015 3.9 
Broad 32 32 xo x 12.6 12.6 99.98740 3.8 


We present results in each of the three standard science bands considered (soft, 
soft+medium and broad), for two values of the significance threshold parameter, 
sigthresh, where sı and s» are significance thresholds of 1 x 10-° and 5 x 1077. From 
the documentation, the significance threshold should be set to s = 1/n,x, which typi- 
cally produces one false detection. The ACIS-S S3 CCD has 1024x1024 pixels, so sı 
should be a sufficiently low threshold. However, to align with the world coordinate system 
(WCS), the CCD image is rotated relative to the bounding box of the image, meaning 
that the image on which wavdetect was run has 1414x1401 pixels, hence the value of 
sy = 1/(1414 x 1401) ~ 5 x 1077. The sky density is computed as the number of sources, 
divided by the field of view, which is estimated as (1024? x 0.49207), where the second num- 
ber is the pixel size in arcsec. Interpreting the sky density as the probability of chance 
alignment, we provide the confidence, CLatign, with which we can reject the hypothesis 


that our source counts could originate from a background source. 
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Extended Data Table 4 | Confidence interval on Chandra ACIS-S count rate. 


Band Energy Count rate (ap.~! ks~') 
68% 90% 
Source Background 
(keV) low high low high 
Soft 0.5-1.2 0.038 0.0010 0.021 0.059 0.013 0.078 
Soft+Medium 0.5-2.0 0.047 0.0021 0.027 0.070 0.018 0.090 
Broad 0.5-7.0 0.047 0.0057 0.024 0.066 0.015 0.086 


The 68% and 90% confidence intervals for the source counts (/V) given in Extended Data 
Table|2} which have been converted to a count rate using the 106.33 ks total exposure time 
for the three standard Chandra energy bands used in this study. Confidence intervals are 
calculated using the method of ref25) which is a Bayesian approach to Poisson statistics 
in the presence of a background that uses a simple prior on the number of source counts 
not being negative. All count rate values are given in units of counts per ks, per 1 arcsec 


aperture. 
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Extended Data Table 5 | Plasma temperature from spectral modelling using XSPEC. 


Plasma Temperature (keV) 
Abundance Model 68% 90% 


Best-fit , , 
low high low high 


Isothermal 0.273 0.198 0.401 0.165 0.653 
Cooling flow 0.317 0.226 0.581 0.187 0.922 


Sola 


Isothermal 0.473 0.231 0.653 0.172 0.772 
Cooling flow 0.594 0.273 0.882 0.193 1.107 


Bulk Earth22 


Isothermal 0.491 0.283 0.661 0.185 0.778 


Photospheric2# 
Cooling flow 0.612 0.330 0.889 0.210 1.111 


The 68% and 90% confidence intervals for the derived plasma temperature are based 
on fits to the full ACIS-S spectrum (plotted in Extended Data Figure |4). Best fits and 
confidence intervals were determined using the C-statistic2 and with no background 
subtraction. For the cooling flow model these are the upper temperature of a range of 
temperatures extending down to 0.08 keV, with the emission measure at each temperature 


weighted by the inverse of the emissivity. 


50 


Extended Data Table 6 | Best-fit X-ray flux from spectral models. 


Flux (10~!° erg s~' cm~?) 


Energy 
68% 90% 
Best-fit , , 
(keV) low high low high 
0.5-7.0 1.78 1.40 2.53 1.20 4.53 
0.3-7.0 1.97 1.49 3.52 1.26 9.40 


0.0136-100 4.34 262 19.1 2.09 218.0 


The best-fit fluxes are computed across three spectral energy ranges using the cooling 
flow model with photospheric abundances?! The first two energy ranges are within the 
wavelength range of our observations and provide a robust lower limit on the flux. The 
third, wider energy band allows for more flux to be carried at higher and lower energies, 
thus providing an upper estimate of the X-ray flux based on our observations. The 68% 
and 90% confidence intervals for the derived flux based on those for the plasma temper- 


ature in Extended Data Extended Data Table|5} 
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Extended Data Table 7 | X-ray luminosity of G29-38. 


X-ray Luminosity (10% erg s~') 


Energy 
68% 90% 
Best-fit 
(keV) low high low high 
0.5-7.0 6.54 5.15 9.30 4.41 16.7 
0.3-7.0 7.24 5.48 12.9 4.63 34.5 


0.0136-100 16.0 9.63 70.2 7.68 801.3 


The best-fit X-ray luminosity (Lx) derived using Lx = 4rd?Fy for the fluxes (Fx) given in 


Extended Data Tablel6] 
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Extended Data Table 8 | X-ray accretion rate of G29-38. 


X-ray Accretion Rate (10° gs~*) 


Energy 
68% 90% 
Best-fit , l 
(keV) low high low high 
0.5-7.0 1.48 1.16 2.09 1.00 3.76 
0.3-7.0 1.63 1.23 2.92 1.05 7.80 


0.0136-100 3.60 2.17 15.8 1.73 181 


The best-fit X-ray accretion rate (Mx) derived using Mx = 2Rwp/(GMwp) with the lumi- 


nosities given in Extended Data Table/7| 
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